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For a capacitor made of a semiconducting carbon nanotube (CNT) suspended above a metallic 
gate, Coulomb correlations between individual electrons can lead to a capacitance that is much larger 
than the geometric capacitance. We argue that when the average spacing n -1 between electrons 
within the low density 1-dimensional electron gas (1DEG) in the CNT is larger than the physical 
separation d between the CNT and the gate, the enhancement of capacitance is expected to be 
big. A recent experiment]- 1 however, has observed no obvious increase of capacitance even at very 
low electron density. We show that this smaller capacitance can be understood as the result of the 
confining potential produced by the potential difference between the source/drain electrodes and 
the gate, which compresses the 1DEG when the electron number decreases. We suggest that by 
profiling the potential with the help of multiple split gates, one can return to the case of a uniform 
1DEG with anomalously large capacitance. 


I. INTRODUCTION 

One of the great features of low-dimensional materi¬ 
als is that their properties can be continuously altered 
using electrostatic gating. Such gating allows one to 
tune the bulk electron density widely and reversibly, 
and thereby probe different regimes of electronic behav¬ 
ior simply by turning an experimental knob. What’s 
more, the material’s response to the applied voltage gives 
information about the electronic compressibility, which 
is in turn a reflection of the nature of the interactions 
and correlations between electrons. For example, for 
a gated two-dimensional electron gas (2DEG), the two- 
dimensional electron density ri2 can be tuned by changing 
the gate voltage V, and the resulting differential capac¬ 
itance C oc dri2/dV of the gate-2DEG system can be 
used as a probe of quantum and many-body correlated 
electron physics. The general promise of capacitance as 
a probe of low-dimensional materials has been appreci¬ 
ated for several decades^®, and has recently been used 
as the basis for a number of dramatic observations in 

graphenS^®. 

One particularly striking observation is that electron 
correlations can lead to an enhancement of the capac¬ 
itance above the normal “geometric” value C g . For ex¬ 
ample, for a low-density 2DEG, the capacitance has been 
predicted to follow the relation^ 

{C/A )- 1 = (< C g /A r 1 - 1.5 n ~ 1/2 (1) 

(in Gaussian units), where A is the capacitor area. The 
second term on the right hand side of Eq. 0 arises as 
a consequence of spatial correlations between electrons, 
which lower the total energy of the electron gas rela¬ 
tive to the uniform state and thereby drive up the ca¬ 
pacitance. While a similar result is often explained in 
terms of the exchange interaction and referred to as a 
“quantum capacitance”, its origin can be understood at 
the level of a classical electron crystal (a Wigner crys¬ 
tal), where positional correlations between electrons re¬ 


duce the total electrostatic energy^ 11 . An analogous re¬ 
sult to Eq. ([I]) has been observed already for a num¬ 
ber of two-dimensional systems, including GaAiP® ; , the 
LaA103/SrTi03 interface®, and graphene in a magnetic 
field™. 

Recently, however, it was pointed out that Eq. ([I]) 
cannot retain its validity at arbitrarily small electron 
density 11 . Indeed, at sufficiently small density that the 
typical distance between electrons n 2 ' becomes much 
longer than the distance d to the metal gate, Eq. ([!]) 
gives a negative result for C. Such negative values of 
the capacitance are forbidden by thermodynamic sta¬ 
bility criteria® (assuming that the capacitance is mea¬ 
sured in the low frequency limit). This failure of Eq. 
0 points to the appearance of new physics in the very 
low-density regime. This physics was outlined for two- 
dimensional systems in Ref. Illl where it was shown that 
at ri2d 2 <C 1 the interactions between electrons become 
strongly screened by image charges in the metal gate, so 
that the electron-electron repulsion takes the form of a 
dipole-dipole interaction. Such screening allows for an 
enormous enhancement of the capacitance over the geo¬ 
metric value in the low-density limit. It also implies a 
re-entrance of the Fermi liquid phase at small ri2, and an 
ultimate saturation of the capacitance at C/A ~ me 2 /h 2 
in the limit ri2 —> 0, where m is the effective electron 
mass®, e is the electron charge, and h is the reduced 
Planck constant. 

So far, however, experiments have been largely unable 
to probe this low-density regime of “dipolar electrons”. 
While some efforts have come closeP, experiments on 
two-dimensional electron systems have been largely hin¬ 
dered by the ability to make very thin devices. The pres¬ 
ence of disorder is also a ubiquitous confounding factor, 
since disorder produces strong spatial modulation of the 
electron density when the electron number is low. 

Recently, however, developments in ultra-clean, gated 
carbon nanotube devices have shown the promise to over¬ 
come both of these difficulties. In such one-dimensional 
devices, a single semiconducting CNT is suspended above 
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a gate electrode, as illustrated in Fig. |TJ and its capaci¬ 
tance is measured at low temperature. Importantly, ex¬ 
isting devices exhibit a very large ratio of the length L 
to the gate distance d, and can have their total electron 
number N tuned incrementally from zero. These condi¬ 
tions together raise the possibility of convincingly prob¬ 
ing the dipole-interacting regime of electron behavior, in 
which the inter-electron separation is much larger than 
d. 


the previous sections, thereby justifying a simple clas¬ 
sical description. We conclude in Sec. [V] by discussing 
the outlook for future experiments. We also discuss how 
the case of a spatially-uniform system can be realized us¬ 
ing a setup with multiple gates, which would allow for a 
significantly larger capacitance. 

II. SPATIALLY-UNIFORM CASE 



FIG. 1. (Color online) Experimental set-up of the CNT de¬ 
vice. A CNT (purple) is suspended above a gate electrode 
(blue) between source and drain electrodes (orange). L is 
the length of the CNT, w is its radius, d is the separation 
between the CNT and the metallic gate, and n~ 1 is the aver¬ 
age distance between electrons. Vds denotes the source-drain 
voltage, which is much smaller than the the gating voltage 
V g . • symbols represent electrons and o’s denote the induced 
positive charges on the gate. 


The one-dimensional nature of a CNT, however, pro¬ 
duces qualitative changes to the capacitance at low elec¬ 
tron density that have so far not been studied theoret¬ 
ically. In this paper we address this problem, describ¬ 
ing the nature of electron correlations across all different 
regimes of the total electron number N and deriving an¬ 
alytical expressions for the differential capacitance. Our 
analytical results are confirmed with simple numeric cal¬ 
culations. We focus everywhere on the case of a semi¬ 
conducting CNT, where electrons have a finite effective 
mass at low density. The case of a metallic CNT, where 
electrons have a linear, Dirac-like dispersion relation, re¬ 
quires separate consideration. 

In describing the behavior of the system, we focus on 
two cases for the experimental environment. In the first 
case (Sec. 0, we imagine that the system is spatially 
uniform, so that electrons are not subjected to any ex¬ 
ternal potential, as in previous work 11 ! However, we find 
a big discrepancy between the predicted results for this 
case and the experimental data in Ref. [I. A considera¬ 
tion of the experimental setup leads us to the second case 
(Sec. Ill), in which electrons are subjected to a confin¬ 
ing electrostatic potential, as might arise naturally from 
the applied potential difference between source/drain and 
gate electrodes. As shown in Sec. m this confining po¬ 
tential is crucial for understanding previously-published 
experimental data. In Sec. [WJ we check our calcula¬ 
tions using a variational method that properly accounts 
for the quantum kinetic energy ignored in Sec. [TT] and 


III The results are found to closely align with those of 


In this section we consider a model in which the CNT is 
not subjected to any external potential, and its electron 
charge is simplified as a ID uniform electron gas which 
forms a capacitor with the metallic gate placed below it. 
The CNT’s radius w is ~ 1 nm, much smaller than its 
length L (which in Ref. CD is 880 nm) or the gate distance 
d (130 nm in Ref. [T]). The linear density of electrons 
along the nanotube is n = N/L , where N is the total 
electron number. The typical distance between electrons 
in experiment is n _1 = 50 — 880 nm. The relevant length 
scale for describing quantum mechanical effects is the 
effective Bohr radius, a B = h 2 /me 2 ss 8.9 nm. This 
large Bohr radius arises in experiment from stretching or 
curvature effects in the nominally metallic CNT, which 
open up a small band gap A = 34 meV, and then the 
effective mass of electrons is small: m ~ mo/170 (where 
to o is the bare electron mass). 

Our general approach to calculating the capacitance 
(as outlined, for example, in Ref. ED is to first evalu¬ 
ate the total energy E of the CNT-gate system and then 
calculate its second derivative with respect to the total 
electron number. In particular, the differential capaci¬ 
tance is given by the formula 


C = 


d 2 E 
e 2 dN 2 


-l 



( 2 ) 


where /i = dE/dN is the chemical potential, which in¬ 
cludes interactions both between electrons within the 
CNT and between electrons and the positive charge in 
the gate. Equation ([2]) implies that capacitance can be 
thought of as equivalent to electron compressibility: large 
capacitance arises in situations where electrons are easily 
compressible. We concentrate throughout this paper on 
the limit of zero temperature, which provides an accurate 
estimate for the capacitance whenever the thermal energy 
ksT is much smaller than the typical interaction energy 
per electron. More generally, for finite temperature sys¬ 
tems one can calculate the capacitance by replacing the 
energy E in Eq. § with the Helmholtz free energy. 


A. High Density Regime: n 1 /a* B 

For an interacting electron system, when the average 
distance between electrons is much smaller than the ef¬ 
fective Bohr radius a * B , the electrons can be considered to 
be weakly interacting and are spread uniformly over the 
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nanotube. Thus, at sufficiently high density, the CNT 
can be modeled as a metallic wire, and the dominant 
contribution to the total energy E comes from the elec¬ 
trostatic energy of the wire-gate system. Such a system 
has total energy e 2 N 2 ln(2 d/w)/L, which corresponds to 
a chemical potential /i = dE/d,N = 2e 2 n \n(2d/w) p so 
that the capacitance is 

c = c » = 2 mSaT (3) 

We refer to this quantity as the geometric capacitance. 

One can note, however, that Eq. ^ gives an ex¬ 
pression for the capacitance that is exact only when 
n 1/y/a* B w. At densities 1 /a* B -C n <C 1 /\/ a B w, 
though the electron gas is still in the weakly interact¬ 
ing regime, the argument inside the logarithm is mod¬ 
ified due to short-range correlations between electrons. 
Such correlations emerge at small enough electron sep¬ 
aration r that the Coulomb repulsion e 2 /r is larger in 
magnitude than the typical kinetic energy ~ h 2 n 2 /m , 
and they produce a truncation of the electron-electron 
interaction at short distances. As a consequence, in this 
density range the argument of the logarithm is replaced 
by dn 2 a* B . Since this replacement produces only a small 
correction, one can still say that Eq. ^ is correct to 
within a logarithmic factor whenever na* B 1. 


B. Intermediate Regime: 1 /a* B n 1/d 


When n <C l/a|j, the typical Coulomb energy ~ e 2 n 
dominates over the kinetic energy ~ h 2 n 2 /m , and the 
electrons assume a correlated state reminiscent of a 
Wigner crystal. (Indeed, such a Wigner-crystal-like state 
has already been reported experimentally in semicon¬ 
ducting CNTs 1 ^.) One can therefore estimate the total 
energy of the electron system by modeling it as a line of 
classical point charges with the gate-screened interaction 
V(r) = e 2 /r — e 2 /\Jr 2 + (2d) 2 . This approach gives 


E = nL 




V(i/n ) 2 + (2d) 2 


( 4 ) 


where the first term inside the braces represents the at¬ 
traction of a given electron to its own image charge, and 
the sum describes the repulsion between electrons. This 
classical description, which ignores the quantum kinetic 
energy of electrons, is justified in detail in Sec. m 
Evaluating Ecp Q in the limit nd 2 1 and taking the 
derivative with respect to n gives 


/Li ~ e 


2 n ln(2 dn) + 0.768n 


I 

4d 


( 5 ) 


so that the capacitance 


C 


L 

2[ln(2 dn) + 1.384]' 


( 6 ) 


This expression implies a marginal, logarithmic increase 
in the capacitance with decreasing electron density, 
C/C g ~ 1 + [ln(l /nw) — 1.384]/[ln(2nd) + 1.384]. Such a 
small correction to the geometric capacitance is similar to 
the situation in 2D, where electron correlations provide 
a small, positive correction to t he g eometric capacitance 
at intermediate electron density 2 11 . This enhancement 
of the capacitance is a manifestation of electron correla¬ 
tions, which allow the electrons to lower their energy for 
a given concentration by avoiding each other spatially, 
and thereby achieve a higher compressibility. 


C. Low Density Regime: n <C 1/d 

In the regime of very low electron density, n<l / d 1 the 
distance between neighboring electrons is much longer 
than the distance between an electron and its image 
charge in the metal gate. As a result, the electron- 
electron interaction takes the form of a dipole-dipole po¬ 
tential: V(r) ~ 2 e 2 d 2 /r 3 . Consequently, the total energy 
can be calculated as 



Therefore, by Eq. ([2]), 

e 2 

—+9.62 e 2 d 2 n 3 (7) 


and 


0.035L 

6 ~ 


( 8 ) 


Thus, when the electron density is made sufficiently 
low that nd 1, the capacitance grows very quickly 
with decreased n, as in the 2D cas^M This large capaci¬ 
tance can be seen as a consequence of the screening of the 
electron-electron interactions by the metal gate, which 
renders them effectively short-ranged. This truncation 
of the interaction, together with the positional correla¬ 
tions between electrons, allows the electrons to be highly 
compressible, and drives the capacitance far above the 
geometric value given by Eq. ([3]). 

Finally, if the electron density is reduced so far that 
n <C a B /d 2 , then the dipolar interaction between elec¬ 
trons becomes too weak to maintain a strong degree of 
spatial correlation between electrons. In other words, 
one can think that the zero-point fluctuations of an elec¬ 
tron in the confining potential created by its neighbors 
become comparable in amplitude to the spacing between 
electrons. In this case, the Wigner-crystal-like correla¬ 
tions between electrons are melted, and the system can 
again be described as a weakly-interacting liquid. In such 
a situation the total energy of the electron system is dom¬ 
inated by the kinetic energy, and the capacitance is pro¬ 
portional to L/na* B . This extreme low density regime is 
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never reached in the experiments of Ref. [T] and is ne¬ 
glected for the remainder of this paper in order to relate 
more closely to experiment. 

The behavior of the capacitance across different elec¬ 
tron density regimes is shown in Fig. [2] The capacitance 
C shows a very weak dependence on the electron density 
when nd 1, while at nd -C 1 the capacitance rises 
quickly with decreasing electron concentration. 



FIG. 2. Capacitance per unit length, C/L , for a CNT with¬ 
out a confining potential, plotted as a function of the electron 
density n multiplied by the gate distance d. The solid line 
shows a numerical evaluation of C using Eqs. @ and CD- 
Dashed lines illustrate the asymptotic behaviors of C in dif¬ 
ferent density limits. At small n, the capacitance is described 
by Eq. |8|. In the intermediate regime, the behavior is de¬ 
scribed by Eq. (|6|. At much larger n, C saturates at the 
geometric value, Eq. Cll- 


Interestingly, the first experimental measurement of 
the capacitance of a single CNT capacitor observed a 
similar trend of increasing capacitance with decreasing 
electron densitjC®. This experiment, however, reported 
only a small increase in the capacitance relative to the un¬ 
correlated state, rather than the large effect depicted in 
Fig. H This absence of large capacitance enhancement 
in this experiment can likely be attributed to the rela¬ 
tively large experimental temperature (T ss 77 K), which 
destroys electron correlations when the electron density 
is low and the typical electron-electron interaction en¬ 
ergy is weak. In particular, for the experiment of Ref. 
HU the temperature was large enough that ksT ~ e 2 /d, 
which implies that electron correlations are lost at all 
n < 1/d and the theory outlined in this section no longer 
holds. Since the increase in capacitance is relatively weak 
at nd > 1 (as shown in Fig. [2]), one can conclude that 
only a relatively small enhancing effect of the capacitance 
should be expected. Disorder effects may also contribute 
to the observed absence of large capacitance at low den¬ 
sity. 

Fortunately, more recent experiments- have performed 
capacitance measurements at much lower temperature 
(4 K), for which strong electron correlations can persist 


down to low densities n <C 1/d. These experiments are 
considered in detail in the following section. 

III. CAPACITANCE IN THE PRESENCE OF A 
CONFINING POTENTIAL 

In the previous section, we derived the dependence of 
the capacitance on the electron density for a CNT system 
without a confining potential, and we showed that C rises 
sharply at small n <C 1/d. A recent experiment how¬ 
ever, showed no such enhancement of the capacitance, 
even when ( N/L)d is made as small as ss 0.15. In this 
experiment, the capacitance was studied by examining 
the Coulomb blockade in the source-drain current (rather 
than by measuring the impedance produced by a small 
AC modulation of the gate voltage, as is typical for 2D 
systems). The authors of Ref. |Tj found that the spac¬ 
ing in the source-drain voltage Vds between successive 
Coulomb blockade conductance peaks was nearly con¬ 
stant, A V g = 31.5 ± 1.5 mV, which suggests a constant 
capacitance C m = e/AV g ~ 46 mn. This spacing was 
consistent over the reported range of electron number, 
1 < N < 17, which corresponds to an electron density 
0.001 nm -1 < n < 0.019 nm -1 . 

Inserting the experimental parameters from Ref. [T| into 
the theoretical prediction of the previous section gives the 
result shown in Fig. [3] As can be seen from this figure, 
this theory predicts a capacitance that is larger than the 
experimentally observed value by as much as 20 times. 
The experimental results also do not exhibit the sharp 
upturn in capacitance at small electron density that is 
characteristic of the theory in the previous section. 



FIG. 3. Prediction for the capacitance C as a function of the 
electron density n using the experimental parameters of Ref. 
IH The solid line is the prediction of the theory for a uniform 
gas (see Sec. E . The dashed line shows the actual capacitance 
Cm = 45.7nm measured experimentally. 


In order to resolve this discrepancy, we consider the 
effects of the confining potential produced by the elec¬ 
trodes. In the experiment, a potential difference V g ss 
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0.5 — 1V is applied between the gate electrode and the 
CNT in order to bring the Fermi level close to the bottom 
of the conduction band. This potential produces signifi¬ 
cant electric fields, as depicted in Fig. [4j which imply that 
the electric potential along the CNT is non-uniform. 

In order to estimate the strength of this confining po¬ 
tential, we solve for the electrostatic potential of a system 
of three coplanar metallic strips held at different volt¬ 
ages (due to the small size of the CNT, we can ignore its 
role in modifying the potential). An analytic solution for 
the potential is presented in Appendix [A] These metallic 
strips are assumed to have infinite length in the direction 
perpendicular to the CNT, but finite width, which in ex¬ 
periment is s ~ 300 nm for the source/drain electrodes 
and L = 880 nm for the gate. This assumption of infi¬ 
nite length is justified by the large ratio of the electrodes’ 
length to their width, which in Ref. [H is > 100. The as¬ 
sumption of coplanarity is justified by the small value of 
d/L. 

If one defines the coordinate x as the position along 
the CNT axis relative to the center of the CNT, then the 
electric potential V has a parabolic maximum at x = 0, 
so that the confining potential can be written 

p 2 

U(x) = -eV(x, d) ~ £7(0) + -^x 2 

where U (0) is a constant and D is a length scale that de¬ 
scribes the strength of the confinement. The expression 
for D is given in Appendix [Aj and for the experiments of 
Ref. dj D « 72- 90 nm. 

When a small number of electrons are added to the 
CNT, the confining potential pushes these electrons to 
occupy a region near the center of the gate electrode, so 
that the electron-populated region has an effective length 
2 L e that is shorter than the total CNT length L. As N 
is reduced, electrons are more easily compressed by the 
confining potential and L e shrinks. This reduction of the 
effective length of the capacitor tempers the increase of 
C at low density that was derived in the previous section. 

In the remainder of this section we calculate the capac¬ 
itance in different regimes of the total electron number, 
taking into account the compression effect produced by 
the confining potential. We focus our discussion here on 
the case where the local electron density is everywhere 
smaller than 1 /a* B and larger than a^/d 2 , which is the 
relevant case for existing experiments. Within this range 
of density, we focus on only two regimes: N/L e d -1 
and N/L e or equivalently (as shown below), N 

( D/d ) 3/<2 and N -C (D/d) 3 ^ 2 . These two regimes corre¬ 
spond, respectively, to the regime of a Wigner crystal-like 
state with unscreened interactions, and the regime of a 
“dipole crystal” with strongly-screened interactions. 


A. High Density Regime: N ( D/d ) 3/ * 2 

Because of the confining potential, the electron den¬ 
sity n varies with position x in order to maintain a 


(a) 




FIG. 4. (Color online) Confining potential acting on electrons 
in the CNT formed by fields between source/drain electrodes 
and the gate, (a) Schematic of the electric field produced by 
the potential difference between the source/drain electrode 
and the gate. The source-drain voltage Vds is neglected here, 
(b) This field produces a confining potential U(x) along the 
CNT which leads to non-uniform electron density, so that the 
electrochemical potential fio = U(x) + fi(n(x )) is constant. 
The effective length 2 L e is defined by n(±L e ) = 0. 


constant value of the electrochemical potential o = 
U(x) + n (n(x)) (see Fig. [4|. In the regime of relatively 
high electron density, electron interactions are essentially 
unscreened by the gate, and the contribution to the chem¬ 
ical potential from electron interactions is given by Eq. 
Q, jn(n(x)) w e 2 n[2ln(2dn) + 0.768]. (Here, the n- 
independent term in /i is ignored, since this just redefines 
the reference point for fi). Thus, the electrochemical po¬ 
tential can be written 

e 2 

/io = -jy^x 2 + n(x)e 2 [2\n(2dn(x)) + 0.768]. 

Rearranging this expression for n(x) gives 

= /i 0 /e 2 - x 2 /D 3 

2 ln(2dn(x)) + 0.768 

A solution can be obtained for n(x) by taking suc¬ 
cessive approximations for the slowly-varying logarithm 
ln(2 dn{x))\ we treat this factor first as a constant of or¬ 
der unity and later replace n(x) inside the logarithm with 
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the calculated average value N/(2L e ). Since 


TV = 2 f n(x)dx , 

Jo 


this approach gives 


L e ~D 



1.237V 

(. D/d ) 3/ 2 


Mo - 


I? 



1.237V 

{D/d) 3 / 2 


1/3 

2/3 


Thus, as announced earlier, N/L e ^> d 1 is equivalent to 
TV {D/d) 3 ^ 2 , and the resulting capacitance 




1.5 


In 


1.231V 

(D/d)3/ 2 


r 


D\ 


3/2 


TV 

{D/d) 3 / 2 _ 


1/3 


( 9 ) 


One can notice that Eq. Q represents a qualitatively 
different dependence of the capacitance on electron num¬ 
ber than the corresponding expression for the uniform 
gas, Eq. ©>■ Rather than weakly increasing as the num¬ 
ber of electrons in the CNT is reduced, Eq. implies 
that the capacitance decreases with decreased TV. This 
decline of the capacitance is a result of the confining 
potential, which causes the effective capacitor length to 
shrink, and this shrinking of the effective length dom¬ 
inates over the weak effect of electron-electron correla¬ 
tions. 


B. Low Density Regime: TV {D/d) 3 / 2 


On the other hand, in the regime where electron- 
electron interactions are strongly screened by the metal 
gate, the effect of electron-electron correlations is large, 
and competes with the shrinking effective length to de¬ 
termine the overall capacitance. In this regime, the con¬ 
tribution to the electrochemical potential arising from 
electron-electron interactions is given by Eq. 0 , so that 
the electrochemical potential is 

e 2 

Mo = + 9.62 e 2 d 2 n 3 {x). 

D 3 


Solving this expression for the electron density gives 


n{x) 


/ Mo-eVV D 3 \ 1/3 

v 9.62 e 2 d 2 ) 


( 10 ) 


Using the expressions TV = 2 n(x)dx and n{L e ) = 0 
gives 


L e = 1.15d 2/5 Z? 3/5 7V 3/5 , 

Mo = 1.32e 2 (T 4/5 T?- 9/5 TV 6/5 , 



FIG. 5. Graph of capacitance C as a function of electron num¬ 
ber TV in the presence of a confining potential. C is expressed 
in units of d{D/d) 3 / 2 and TV is normalized by {D/d) 3 / 2 . The 
solid line shows the numerical result for capacitance. The 
dashed lines represent the asymptotic behaviors described by 
Eqs. Q and (111. 


so that, as announced earlier, the condition N/L e -C d 1 
is equivalent to TV <C ( D/d) 3 ^ 2 . The corresponding ca¬ 
pacitance C = e 2 {d^ 0 /dN )~ 1 is given by 


C 


0.629d 



3/2 


TV 

{D/d) 3 / 2 _ 


-1/5 


( 11 ) 


Thus, in the low density regime, the system exhibits 
an increase in the capacitance at small electron density, 
as in the case of a uniform electron system, due to the 
high electronic compressibility associated with the strong 
screening of the electron-electron interaction. However, 
the growth of the capacitance is largely muted by the 
shrinking effective length, and retains only a very mild 
dependence, C ~ TV -1 / 5 . 

A plot of C'(TV) is presented in Fig. [5] showing the dif¬ 
ferent regimes of behavior. Notably, the capacitance no 
longer varies monotonically with TV, first decreasing with 
TV in the low density limit and then increasing again at 
high density. 

This theoretical result is compared to the experimen¬ 
tal measurements of Ref. [T] in Fig. [6] The results are 
qualitatively consistent with each other, suggesting that 
the confining potential provides an adequate explanation 
of the small, relatively constant capacitance seen in ex¬ 
periment. A more accurate comparison to experiment 
may require more careful modeling of the electrostatic 
potential created by the experimental setup. 


IV. CORRECTIONS TO CAPACITANCE DUE 
TO QUANTUM EFFECTS 

So far, our calculations of the capacitance have em¬ 
ployed a purely classical approximation, in which the 
quantum kinetic energy of electrons is ignored. This 
approximation is generally justified in the strongly- 
interacting regime, where quantum effects provide only 
a small correction to the total energy. Nonetheless, in 
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FIG. 6. A comparison between the predicted behavior of 
C(N) and the experiments of Ref. [L The solid line is the 
prediction of the theory derived in Sec. m The dashed line 
represents the experimentally-measured capacitance C m . 


this section we carefully justify our previous results by 
performing a calculation in which quantum mechanical 
effects are properly taken into account. 

Our goal is to accurately estimate the ground state 
energy E of the interacting electron system as a func¬ 
tion of the total electron number N ; the capacitance can 
then be found by Eq. © In other words, we seek an 
accurate estimate of the lowest-energy eigenvalue of the 
Hamiltonian 


* = -£ 

i 


2 m 




(The additional interaction between an electron and its 
own image charge does not enter the capacitance, and 
can be ignored.) In order arrive at an estimate of E, we 
use a variational method, in which the expectation value 
of the Hamiltonian, (H)\ = is calculated for 

a wide set of variational wavefunctions {|d>^)}, where A 
denotes some variational parameter that labels the dif¬ 
ferent wavefunctions in the set. The ground state energy 
E is associated with the minimum value of (H)\ across 
all possible values of A. 

Following Ref. 1201 we use for our variational wavefunc¬ 
tions the eigenstates of the exactly-solvable Calogero- 
Sutherland model (CSM)p2 so that the variational pa¬ 
rameter A describes the CSM interaction constant; A = 1 
corresponds to the non-interacting limit, while A —> oo 
corresponds to a perfectly-ordered crystal. In the limit 
of large system size, the CSM eigenstates have a kinetic 
energy per particle given by" 


£fc(n, A) 


7 r 2 h 2 n 2 A 2 
6m 2A — 1 ’ 


while the interaction energy per particle is 


Gni(^b A) — rj 


V (r)g( A, r)dr. 


Here, g(X,r) is the pair distribution function associated 
with the eigenstate While g(X,r) has been studied 


in some detail analyticall y — I— ( closed-form expressions 
are known only for the special values A = 1/2, 1, 2, and 
oojm Thus, following Ref. [201 we evaluate e int (n,X) at 
these special values, and then interpolate to other values 
of A by making a third-order polynomial fit as a function 
of 1/A. 

The ground state energy per electron, e(n), for a given 
density n is equated with the minimum of the function 
ejfe(n, A) + €int{n,X) across all values of A. The chem¬ 
ical potential is then calculated as fi = d(ne(n))/dn. 
For the case without a confining potential, the capac¬ 
itance is equated with e 2 L(dn/d/r), as implied by Eq. 
([2]). For the case with a confining potential, the capac¬ 
itance is found by calculating the electrochemical po¬ 
tential associated with a given total electron num ber, 
Ho = /i(n(x)) + e 2 D 2 /x 3 , as explained in Sec. |III[ and 
then setting C = e 2 (dN/dfj,o)- 

The resulting capacitance is plotted in Fig. [7] for the 
range of experimental parameters explored in Ref. [lj As 
can be seen, the results closely follow those derived from 
the simple classical description of Sec. [TT] and m sug¬ 
gesting that quantum effects indeed provide only a small 
correction within the experimental range of electron den¬ 
sity. One can also see in Fig. [8] the non-monotonic depen¬ 
dence of the correlation strength on the electron density. 
Indeed, as the electron density is first reduced beginning 
at large density, the electrons become more strongly cor¬ 
related as the typical ratio e int /efc ~ (e 2 n)/(H 2 n 2 /m) ~ 
is increased, and correspondingly the value of 
A associated with the variational wavefunction increases. 
On the other hand, at small enough density that n <C 
a* B /d 2 , the interactions between electrons is dipolar, so 
that the ratio e,„t/efc ~ (e 2 d 2 n 3 ) / (h 2 n 2 /m) ~ nd 2 /a* B . 
Consequently, the correlations between electrons become 
weak again as the density is reduced further, and the 
value of A falls. 



FIG. 7. (Color online) Plot of the capacitance C against the 
number of electrons N. The two solid lines of group A are re¬ 
sults for the case without a confining potential. The thick line 
(blue) is the result using the variational method presented in 
Sec. |IV[ while the thin line (red) is calculated using the classi¬ 
cal approach in Sec. [TT] The two solid lines of group B have the 
same meaning but are for the case with a confining potential. 
The dashed line (blue) indicates the experimentally-measured 
value. 
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(a) 



FIG. 8. (Color online) Variational results for the chemical po¬ 
tential p and the variational parameter A as a function of the 
dimensionless concentration na* B at d = 15 a* B . (a) The thick 
solid line (blue) represents /r(n) in units of e 2 /a* B . The thin 
solid line (red) corresponds to the result obtained by the clas¬ 
sical model in Sec. [TT] (b) The variational parameter A associ¬ 
ated with the minimum energy wavefunction at a given den¬ 
sity n. Vertical dotted lines denote different regimes of elec¬ 
tron density; from left to right: the low-density weakly inter¬ 
acting regime at n <S a B /d 2 , the dipole-interacting Wigner- 
crystal-like regime at a B /d 2 -C n <C 1/d, the unscreened 
Wigner-crystal-like regime at 1/d < n < 1 /a B , and the 
high-density weakly interacting regime at l/a* B <C n. In the 
crystal-like regimes, the difference between variational and 
classical results for fi is small and A is relatively big, denoting 
strong correlations. In the weakly-interacting regimes, the 
difference in ^ grows and A approaches unity. 


V. CONCLUSION 


In this paper we have studied the ID correlated elec¬ 
tron gas in a gated carbon nanotube device. We cal¬ 
culated the capacitance of the system in different den¬ 
sity regimes, which exhibit different electronic behav¬ 
iors. For a spatially uniform system, the electrons oc¬ 
cupy a weakly-interacting phase at high electron density, 
and correspondingly the capacitance is constant and rel¬ 
atively small. As the electron density is lowered, the 
electrons acquire strong, Wigner crystal-like correlations, 
and the capacitance becomes weakly enhanced. At suf¬ 
ficiently low electron density that nd <C 1, the electron- 
electron interactions are strongly screened by the metal 
gate, and the capacitance increases sharply. 

These predictions for the uniform system, however, are 
qualitatively different from existing experimental mea¬ 
surements, which show a capacitance that is roughly in¬ 



(b) 



FIG. 9. (Color online) Possible realization of constant po¬ 
tential along the CNT by using multiple split gates set at 
different gate voltages, (a) A schematic of the device setup, 
taken (with permission) from Ref. [I] (b) By adjusting the 
voltages on the multiple gates, one creates a potential 1 / 2 ( 2 ;) 
which compensates the confining potential Ui(x) caused by 
the source/drain-gate fields, resulting in a uniform potential 
V g = Ui(x) + 1 / 2 ( 2 ;) across most of the CNT. 


dependent of electron concentration and much smaller 
in magnitude. This discrepancy can be understood by 
accounting for the effects of a confining electrostatic po¬ 
tential that tends to push electrons toward the center of 
the CNT, and thereby reduce the effective length of the 
capacitor at small electron density. This shrinkage of the 
effective capacitor length greatly reduces the capacitance 
at small electron density, and results in an estimate for 
the capacitance that is much more consistent with exper¬ 
iment. 

If, in the future, the confining potential could be elim¬ 
inated, the dramatic increase in capacitance at low den¬ 
sity would be restored, and this would permit more sen¬ 
sitive probing of correlated “dipolar electron” physics. 
Further studies should thus consider how to realize this 
reduced confining potential in the experimental environ¬ 
ment. One possible way is to use multiple gates with dif¬ 
ferent voltages on each one, as depicted schematically in 
Fig. |9j in order to tune the electrostatic potential along 
the CNT. Such devices have already been fabricated^, 
but careful attempts to realize a constant potential have 
not yet been reported. Such a constant potential will be 
achieved when a parabolic potential maximum simulated 
by multiple gates precisely balances out the parabolic 
minimum created by the source/drain-gate potentials. 

Finally, as noted above, we have so far not consid¬ 
ered the effect of finite temperature, while the experi- 
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ment of Ref. Q]was performed at T = 4 K. Nonetheless, 
throughout the experimental range of electron density, 
the typical electron-electron interaction energy e is such 
that e/fcs > 96 K> 4 K, and so our zero temperature 
approximation is justified. 
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Appendix A: Potential Produced by Three Coplanar 
Metallic Strips 


To calculate the potential V produced by three copla¬ 
nar metallic strips placed next to each other, we first 
consider them as separated from each other by a certain 
distance b — a and then take the limit of b — a = 0 to 
get the final result. The middle strip is set at a potential 
V g while the side ones are at V = 0. The electric field 
is independent of the longitudinal coordinate (perpendic¬ 
ular to the CNT and parallel to the gate electrode), so 
that solving for the potential can be reduced to a two- 
dimensional electrostatic problem, as shown in Fig. |10| 
Therefore along the x axis, we have a field 


E x (x) 


2 p(x')dx' 


where p{x') is the surface charge density on the metallic 
strips at x = x'. 



FIG. 10. A schematic graph of three coplanar metallic strips 
separated by a distance b — a and held at different potentials. 


Then using Kramers-Kronig relations, we can define 
a complex function /(C), C = x + ' lz i which satisfies 
Ref(x) = E x (x) and Im/(x) = —27 t p{x). Since 


p(x) = 0, a < \x\ < b or |x| > c 

E x {x) = 0, b < \x\ < c or |a;| < a 


we guess /(C) 


A 

yj (x 2 —a 2 )(x 2 —b 2 )(x 2 — c 2 ) 


. Then, 





A 

a 2 )(b 2 — x 2 )(c 2 


dx. 


Taking b = a(l + £), and expanding for £ <C 1, we have 


V n 


A 7T 
a\Jc 2 — a 2 8 


A 


8 a ■/ c 2 — a 2 V g 

7T 


Using the uniqueness theorem, we know this is the only 
solution to this configuration. Then on the middle strip, 
the charge density is 


p{x) 


4ay/c 2 —a 2 V g 
7T 2 (a' 2 — x' 2 )y/c 2 —x 2 ’ 


|ar| < a 


and the electric field in the vertical direction, E z (x, 0 + ) = 
2iTp(x), is 


E z (x, 0+) 


8Vc 2 - a 2 V q 8Vc 2 - a 2 V q ( 1 


nac 


nac 


2 c 2 


where i<a. So near the surface of the middle metallic 
strip at a height of d, the electric potential is 

V(x,d) « V g - E z (x,0 + )d = V(0) - x 2 


where 


U(0) 

D 


t r 8dV c 2 — a 2 V g 

9 nac ’ 

nace 2 _1_ 

_8 dVc 2 -a 2 eV g l/a 2 + l/(2c 2 ) 


If we denote the width of the middle strip as L and 
that of the side ones as s, then a = L/ 2, c = L /2 + s, 
and 


U(0) 


= ^ S - 


32 dyj 

7zL(L+2s) j 


D 


n(L/s+2)e 2 _1_ 

MdyJL/s+leVg 2+l/(l+2s/L) 2 


L. 


Inserting the experimental parameters from Ref.jl](c? = 
130 nm, L = 880 nm, Vg « 0.5 — IV, and s « 300 nm) 
gives D ss 72 — 90 nm. 
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